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We obtain a closed-form analytical expression for the zero temperature Fourier transform of the 
2fcF component of the density-density correlation function in a Luttinger liquid with different spin 
and charge velocities. For frequencies near the spin and charge singularities approximate analytical 
forms are given and compared with the exact result. We find power law like singularities leading to 
either divergence or cusps, depending on the values of the Luttinger parameters and compute the 
corresponding exponents. Exact integral expressions and numerical results are given for the finite 
temperature case as well. We show in particular how the temperature rounds the singularities in 
the correlation function. 
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I. INTRODUCTION 

Interacting one-dimensional systems have been proven 
to exhibit exceptionally rich physicsi such as powerlaw 
decay of the correlation functions with non-universal ex- 
ponents depending on the interaction and exotic phe- 
nomena such as spin-charge separation. Some of these 
phenomena have been identified experimentally in the 
various realizations of one dimensional systems such as 
organic conductors^, nanotubes^i^ or quantum wires^. 
The standard paradigm describing the low energy prop- 
erties of these systems is known as Luttinger liquid 
(LL) theoryi^ii While the correlation functions in LL 
theory are most conveniently represented in space-time 
variablesi*^, often experimental results are interpreted 
more directly in terms of the Fourier transform into mo- 
mentum and frequency space. Because of the branch- 
cut singularity structure of the real space zero temper- 
ature correlation functions, obtaining the Fourier space 
representation is not always straightforward and requires 
care in evaluating the integrals. Perhaps this is why ex- 
plicit forms of the momentum and frequency dependent 
response functions have been comparatively slow to enter 
the literaturei2ii£ 

The Fourier transforms of the zero temperature sin- 
gle particle Greens function were computed in the early 
1990siiii^ while their finite temperature versions followed 
a few years later J^ii^ Whereas in these references closed 
analytical forms for the spinless Luttinger liquid are pro- 
vided, the spinful case appears to be much more diffi- 
cult, and besides numerical results, only a few analytical 
results are knowni^. Moreover, even at zero tempera- 
ture, an exact closed- form expression for the 2fci? part 
of the density-density correlation function for the spin- 
ful case appears to still be absent from the literature^S; 
some special cases of a Luther-Emery liquid with two 
gapped modes are studied in Ref. [l6j. In this paper 
we provide an exact expression for the Fourier trans- 
form of the density-density correlation function in the 



realistic case of different spin and charge velocities and 
essentially arbitrary coupling constants in the spin and 
charge sectors. These results have potential implications 
for several experiments in one dimensional systems. Let 
us mention for example Coulomb drag between quantum 
wiresiiii^ and measurements of the voltage noise on a 
metallic gate in close proximity to quantum wire.-iS. In 
the case of strong interactions where the spin incoher- 
ent regime can be obtained j^^i^^ i^ the finite temperature 
Fourier transform of the 2fcF density correlations contain 
important information about the LL to spin-incoherent 
LL crossover. 

This paper is organized as follows. In Sec. |TT] we in- 
troduce the notation and conventions we will use to de- 
scribe the spin and charge sectors (including the correla- 
tion functions) of the LL. In Sec. Illll we obtain an exact, 
closed-form expression for the Fourier transform of the 
zero temperature density-density correlation function in 
the important case of different spin and charge veloc- 
ities. We present approximate zero temperature results 
near the spin and charge singularities in Sec. II VI and finite 
temperature results in Sec. |Vl We summarize our main 
points in Sec. IVII A few useful results and expressions 
are relegated to the appendices. 



II. THE SPINFUL TOMONAGA-LUTTINGER 
MODEL 

The low energy properties of a ID system of spinful 
fermions can be studied with the following Hamiltonian 
in bosonizcd formii 
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/)cr) is a bosonic field representing charge (spin) 
collective mode oscillations, 9^ is the dual field satisfying 
[4>u{x),dxO^{x')] = iiT5^^5(x—x'), Vi, are the propagation 
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velocities of these modes, and K^, their stiffness constants 
(in this paper we take H = 1). The Hamiltonian ([1]) be- 
comes SU(2) invariant for the special value K„ = 1. The 
fact that the charge and spin fields commute leads to 
the well-known effect of spin-charge separation, a conse- 
quence of which is the factorization of certain correlation 
functions into a product of spin and charge components 
when expressed as a function of space and time. It turns 
out that for the line Kp -f = 1 in parameter space 
this factorization also occurs in Fourier space. 

In the bosonic language, the density operator has the 
representation 

p = Po dxfpp + 2/9o cos \^2kFX — \/2(f>pj cos V2(j)a- 

+ 2pQ cos (Akpx - 2%/2(/)p) . (2) 



The first term is the average density po, the gradient term 
represents the density oscillations with zero momentum, 
and the third and four terms are the 2kF and Akp parts of 
the density fluctuations, respectively. The decomposition 
([2]), in turn, leads to an analogous decomposition of the 
Fourier transformed density-density correlation function, 

x(9,w) = xo(g,w) + X2kF{q,uj) +x4kF{q,i^), (3) 

where we have neglected higher order subdominant con- 
tributions. Our main interest is in the computation of 
the 2kp part for arbitrary values of the collective mode 
velocities and of the Luttinger parameters K^. 

In coordinate space and imaginary time, the time or- 
dered density-density correlation function is given by 



{TrP2kf {x, T)p2kF (0, 0)) = pI cos(2fcFa;) J| 



(^) 



'^-p,<^ l^sinh (^1^ {x + iVyT — lasignr)^ sinh ^^3^ {x — iVj^r + iasignr) 
I 
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where (3 is the inverse temperature, and a is a short- 
distance cutoff of the order of the inverse bandwidth.'^" 
We shall follow the standard procedure to find X2fcF {li 
the Fourier transform of Eq.(|3]), and then analytically 
continue the result to get the retarded function. In spite 
of the simple factorized expression ([3]), a closed analytical 
form of its Fourier transform has been obtained only for 
the case where both spin and charge velocities are equal, 
or both Luttinger parameters are equal to onci^ For the 
SU(2) invariant case with arbitrary Kp, the double inte- 
gral resulting from the Fourier transform can be reduced 
to a single integral, but it cannot be evaluated in closed 
form. However, if we set ia = in the denominator of 
(HI), the double integral can be performed at zero tem- 
perature for arbitrary values of Vi, and and K^, inside the 
parameter regime Kp + K„ < 2. Last, these results can 



be extended to the region 2 < Kp + K„ < 4 where the 
Fourier transform bears the same functional form than 
for Kp + K„ < 2 plus a constant that depends on a. In 
particular, the singular behavior near k = Up.o-w is given 
in both regimes by and (fTB|) 

III. EVALUATION OF X2fcp(g,t<^) AT ZERO 
TEMPERATURE 

The Fourier transform of ([4]) at zero temperature can 
be written as 



X2feF (9: = 2 [x2feF (9 + 2fcF, iv) + X2feF (9 ~ 2fcF, iv)] , 

(5) 

with 



X2k^{Q^^^) -Po d.^dT —, (6) 

where we have taken ia — > in the denominators^^. To proceed, the denominators must be combined. One way to 
do this is through the introduction of another integral by making use of Feynman's parametrization formula!^ 



A-Bb r(a)r(5)7o [wA+il-w)Bf+''' 
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where F is the gamma function, and we let A = x"^ + v'^t'^, B — + w^r^, a = Kp/2 and b = Ka/2. With this 
formula we find 



X2fcFW>«f^) = / fixdr / dw- 



r 



(^)r(^) 



By interchanging the integration orders, one can evaluate the integrals in x and t as 



(8) 



(^)r(^) 



dw \ 



w 



-Kp + -ff<^-l 



(9) 



where in the last step the restriction Kp + Ka- < 2 holds in order to obtain a finite result. This restriction is an artifact 
of letting ia —^ in the denominators. Had ia been kept finite, we would have obtained a more general result, but 
we would have been unable to perform the integrations so simply. However, it is easy to extend these results to the 
region 2 < Kp + Ka < 4. This can be done by noting that the piece that diverges in the limit a does not depend 
on k and h'. To extract this constant we write the imaginary exponential factor in the integrand of ([5]) as the sum 
[g-^l^aj-i/r) — 1] -|- 1; in the first term in square brackets and since the divergence has been subtracted, one can safely 
take the limit a — > 0, which leads to Eq. ([9]) but now restricted to 2 < Kp + K„ < 4. The second term just adds a 
cutoff dependent constant. Finally, the integral in w can be performed, and after some algebra we find 



4.pg(f)-^+-'F(l-^) 



^ Kp Kp + K„-1 ^ Kp + K, Kp + K, ^ 



(10) 



r 



where Fi is Appell's hypergeometric function of two 
variable o^^i^^ and r = 1 — v^/v^. Interestingly, when 
Kp + Ka = 1 these expressions greatly simplify. The de- 
nominator in the integrand of equation ^ equals unity, 
and the integral reduces to a standard Gauss hypergeo- 
metric function. Thus with Kr, + /^o- = 1 one arrives at 



the simple result 



2 2\ 



(11) 



which shows that for this line in parameter space, the fac- 
torization of the correlation function that represents spin- 
charge separation is also obtained in Fourier space. How- 
ever, the latter factorization is non-trivial in the sense 
that each factor in X2fcF Fourier trans- 

form of the corresponding factor in Eq. (|4]). 



tion function 'x^j^^{q,uj). This, in turn, is related to the 
Matsubara function by means of the analytic continu- 
ation xlk^iq^^) = xlkM^^^ = ^ + iS)- This can be 
readily performed by putting a branch cut for the power 
function on the negative real axis, with the convention 
that it is continuous from above. Consistently, the Ap- 
pell function Fi(a, /?, /?', 7; zi, Z2) has branch cuts in zi 
and Z2 along the interval (l,oo), where it is continuous 
from below 1^ If we assume that Vp > v^, and restrict 
ourselves to w > 0, the limit (5 — > is straightforward, 
except for v^-q < to < Vpq, where the argument Z2 of Fi 
lies over the branch cut. As 5 ^ this line is approached 
from below and therefore, we employ the continuity of the 
function. Thus, we obtain the real frequency version of 

Eq. (nni), 



IV. APPROXIMATE EXPRESSIONS FOR 
X^kF{q,tj) NEAR THE SPIN AND CHARGE 
SINGULARITIES 



Physical properties of the system, i.e., measurable 
quantities, can be extracted from the retarded correla- 
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(12) 



where is the step function, a form that is well suited 
to plotting. 

Many times, useful information can be extracted from 
the behavior near singular points of the correlation func- 
tions. Whereas in one dimension true long-range order 
does not existji for certain classes of models, correlation 
functions decay algebraically for long distances, display- 
ing what has been called quasi long-range order. Alter- 
natively, they show an algebraic singular behavior when 
translated into momentum space. In our result (|12p sin- 
gularities located at w = ifp? and uj = ±Va-q are present, 
and represents charge and spin density fluctuations re- 
spectively, a manifestation of the spin-charge separation 



in a LL. In what follows, we give the explicit behavior of 
p2| near these singular points. 

The Appell function Fi{a, (3, f3' zi, Z2) is defined 
through a double series in zi and Z2 within the conver- 
gence region \zi\ < 1 and I22I < 1-^^ It also has several 
singular points, for instance (21,22) = (0, 1), (0,oo), etc. 
Several transformations permit one to obtain its analytic 
continuation^^ outside this domain by separating Fi into 
two terms: a singular term and a term with a well defined 
expansion around the singular point. This separation al- 
lows us to extract the behavior near the charge and spin 
singularities: 



(I) 



for 



±Vpq 



A, 



\uj — vtq ' 



I D I 2 2 2 

+ Bp\uj - v^q 



-1 -i 

e 



, (13) 



X2fcF («:W) 



Kp+K^ 



for : 



±v„q 



'? 2 I 



- + Kp 



_^_„(i|.+^^_l) B(u--vW) ^ 5^ 12 _ ^ 



2„2| 



\ 9 9 9 \ 

\uj^ — v'lq'^ 



, (14) 



where Apj^ and Bp/o- are given in Appendix [Al 

More schematically, the singular behavior is given by 

X§fe,(<?,^)~ |cc>2-i;y|*+^'-i-f const (15) 

for UJ w i^pQ and u w iva-q respectively. In Fig. [1] we 
show plots of the real and imaginary parts of X2fcF (^i ^) 
(as defined in Eq. ([5])) as a function of w, for a fixed 
value of (7 = and for a ratio of velocities Va/vp = .4. 
The frequency is expressed in units of the charge energy 
Ep = Vpkp. The singular behavior near charge and spin 
singularities (fT5|) and (|14p is also shown. There will be a 
divergence near the charge singularity for -I- K^r < 1 



and a dip otherwise. In cases where a divergence is 
present, it can be accompanied by a change of sign when 
the frequency goes through the singularity (w < Vpq to 
Lu > Vpq) as in shown Fig. [ijb) (real part) [contrary to 
Fig. [IJa)]. This is due to the imaginary exponential fac- 
tors in Eq. (fT5|) . and it is easy to see that the condition 
to be fulfilled for a change of sign is 4- Ka — 1 < — ^ . 
There are analogous results for the spin singularity at 
u! « iva-q- Notice that the presence of prefactors and 
additional constants in Eq. and may change 

dips into peaks of finite height, as is the case of the sin- 
gular behavior near the spin singularity of the real part 
of X2kw in Fig- [Ha), and near the charge singularity of 
the imaginary part of X2fcF in Fig. [Ijc) and (d). 
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FIG. 1: (Color online) Real and imaginary parts of X2k-p {q, 
as a function of frequency (filled red lines) for q = 0, Va/vp = 



.4, and (a) Kp = 1, 



.4, (b) Kp 



.4, K„ 



1 and (d) Kp 



.8, K^ 



.2, (c) 
1. The last 



case corresponds to the repulsive Hubbard model away from 
half filling. The comparison with the approximate expressions 
from Eqs. (|13p and (|14p are also shown (dashed blue lines). 
The charge energy is defined as Ep = kpVp 




FIG. 2: (Color online) Imaginary parts of X2kp{q,^) as a 
function of frequency for q = 0, v^/vp = .1, Kp — .5, K^ = 
1. In the inset we observe the suppression of the correlation 
function under the effect of temperature (we show its behavior 
for fixed values of tj). 



B 2 - 




Lu/Ep 



FIG. 3: (Color online) Imaginary parts of Xzfcp ^) ^ ^ 
function of frequency for q = k-p/2, Va/vp = .2, Kp = .4, 
Ka = .2. The four peak structure is the effect of taking 
9/0. 



the charge singularity is always finite, and the spin sin- 
gularity can be divergent for Kp < 1/2. In particular, 
Eq. (U)) represents the low energy effective Hamiltonian 
for the repulsive Hubbard model away from half-filling, 
in which case the spin parameter flows to the fixed point 
K* ~ 1 and \ < Kp < I. In this situation X2/cf shows 
no divergences in the spin and charge singularities, as it 
is depicted in Fig. [Hd). 



In the figure the four possible cases are presented. In 
Fig. [Tfa) we observe a divergence in the charge singu- 
larity in the imaginary part; in (b) both singularities are 
divergent; in (c) the divergence is in the spin singularity 
and in (d) there are no divergences. One especially im- 
portant case is the SU(2) symmetric line K„ = 1, where 



V. X2fcF ('7.'^) AT FINITE TEMPERATURE 

The finite temperature case is important as it tells how 
the singularities at w = ±W(t9, ±i'p9 will be rounded with 
temperature. Also, if v^jvp <ti 1 the spin-incoherent 
regimeSS, (defined as Wo-fep ^ kgT <^ Vpkp) can be 
approached from temperatures below the spin energy. 
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Vakp. As we will see below, in this limit the 2k-p cor- 
relations (which may dominate the experimental signal 
at zero temperature) are rapidly suppressed with tem- 
perature. The presence or absence of 2fcF density corre- 
lations at a given temperature has imporant implications 
for Coulomb dragi^ and voltage noise^ on a metallic de- 
vice in close proximity to a quantum wire. Moreover, the 
temperature dependence of these quantities can reveal 
spin-incoherent Luttinger liquid behavior 

Unfortunately, we were not able to obtain a simple 
analytical form for the Fourier transform of the finite 
temperature density correlations ([J). However, we were 
able to extract some approximate results and study the 
Fourier transform numerically. 

We study the temperature dependence of the Fourier 
transform using the result jBSj, 

X Imxf (g -q',u;-e') luixfiq',e'), (17) 

wherei 

Xa [Qy^J — 

Va \ Va J 

f .u}~Vak K„\ 
X B —I , 1 

f .uj + Vak Ka 

X B —I , 1 (18) 

V 47rr 4 ' 2 ) ^ ' 

and B is the beta function. A similar expression exists 
for Imxp'"(g,w). 

The behavior of the imaginary part of X2k-p{<lT^) ^ 
function of uj for fixed q is shown in Figs. [5] (for q = 0) 
and [3] (for q ^ 0) for two different sets of parameters. 
Notice that the four peak structure of Fig. [3] is the effect 
of taking g ^ 0. We observe that the main effect of finite 
temperature is 2-fold: i) there is a suppression of the 
whole 2fcF correlation function, and ii) temperature just 
rounds the singularities at w = ±v^q, ±Vpq for Uo-gl < 
nTvp/vcr, while for larger jw— the temperature effects 
are minimal. 



VI. SUMMARY 

We have given an exact closed form expression for the 
zero temperature Fourier transform of the 2kp compo- 
nent of the retarded density-density correlation function 
in a Luttinger liquid with different velocities of spin and 
charge oscillations and arbitrary stiffness constants. Ad- 
ditionally, we have found approximate expressions near 
the collective spin and charge mode singularities that es- 
sentially take the form of power laws, whose exponents 
depend in a simple manner on Kp and K^. We also 
compared these approximations directly with the exact 
result. 

We were not able to find an exact result for the fi- 
nite temperature case, but we were able to evaluate the 
Fourier transform numerically and determine some ap- 
proximate results. One important result of the analysis 
is that the 2k-p oscillations are dramatically (exponen- 
tially) suppressed with temperature when the spin and 
charge velocities are very different. This has implications 
for observable quantities that depend on the 2fcF density 
correlations, such as Coulomb drag or voltage fluctua- 
tions on a metallic gate proximate to a quantum wire. 
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APPENDIX A: COEFFICIENTS OF 
EXPANSIONS NEAR CHARGE AND SPIN 
SINGULARITIES 

The coefficients ^p/o- and i?p/o- appearing in and 
(fT4)) are given by 
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r (i -K,- ^) 



^p = — \ ' ' F\Kp + K„-\, ^'^^" \ Kp + K,-l;r] , (Al) 



r ( 1 - ^"t^" 1 r ( ^ + - 1 



r(i_^^) 



r(^)r(i^p + if,-i) V 2 2 



I 

where F is the Gauss Hypergeometric function (some- where the outer sum can be evaluated with standard com- 
times written 2^1)- plex integration: 



APPENDIX B: GENERAL FORMULA FOR 
FOURIER TRANSFORM OF 2fcF RESPONSE 
FUNCTION 

The 2fcF part of the density-density correlation func- 
tion ([3]) can be expressed as 

x(2:,t) = Xa(a;,r)xp(x,T), (Bl) 

where we use definition of the Fourier transform 

dx \ dTe'^-^-''-\{x,T). (B2) 
-00 Jo 

to write x('7; ^J') as a convolution of terms 



Using the spectral representation 

1 [ ^ Imx''^* 



z — e 



(B3) 
(B4) 



finds 



1 s-^ f dq' f dei f de2 



iv — iuJn — El iijJn — <^2 



, (B5) 



1 ^ 1 1 

f3 ^ iv — iuJn — ei — £2 



"B(e2) - nB(-ei) 



iiy — €2 — ei 



(B6) 



Therefore, one finds 



Xiq,i'^) = - / ^ / 



dq^ f dei f de2 713(62) ^ ?iB(-ei) 
2^ 



TT J TT If — ei — €2 

xlmxr(g-<z',ei)lmxf (^',£2), (B7) 



from which the analytical continuation iv lu + iS can 
readily to be done to yield 



Imx(g,cj) = ^ j J ^[72B{e' ~ uj) ~ nB{e')] 

X ImxTiQ -q\^-e') Imxf (9', e')- (B8) 
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